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Abstract 

We calculate the 0(a 2 ) massive operator matrix elements for the twist— 2 operators, which 
contribute to the heavy flavor Wilson coefficients in unpolarized deeply inelastic scattering 
in the region Q 2 3> m 2 , up to the 0(e) contributions. These terms contribute through 
the renormalization of the O(of) heavy flavor Wilson coefficients of the structure function 
F2(x,Q 2 ). The calculation has been performed using light-cone expansion techniques with- 
out using the integration-by-parts method. We represent the individual Feynman diagrams 
by generalized hypergeometric structures, the e-expansion of which leads to infinite sums 
depending on the Mellin variable N. These sums are finally expressed in terms of nested 
harmonic sums using the general summation techniques implemented in the Sigma package. 



1 Introduction 



The heavy flavor corrections to deeply inelastic scattering constitute an important part of the 
structure functions in the lower x region, cf. [1]. The current world data for the nucleon structure 
functions F 2 ,d (x, Q 2 ) reached the precision of a few per cent over a wide kinematic region. 
Therefore both for the determination of the QCD scale Aq CD and the detailed shapes of the 
partonic distribution functions the analysis at the level of the 0(a^) corrections is required to 
control the theory-errors on the level of the experimental accuracy and below [2]. In a recent 
non-singlet analysis [3] errors for a s (M§) of 0(1.5 %) were obtained extending the analysis 
effectively to N 3 LO. In the flavor singlet case the yet unknown 3-loop heavy flavor Wilson 
coefficients prevent a consistent 3-loop analysis. Due to the large statistics in the lower x region 
one may hope to eventually improve the accuracy of a s (M|) beyond the above value. 

The heavy flavor corrections to F 2 ,d (x, Q 2 ) were calculated to 2-loop order in the whole 
kinematic domain in a semi-analytic way in x-space in Refs. [4]. A fast implementation for 
complex A^-space was given in [5]. In the range of higher values of Q 2 one may calculate the 
heavy flavor Wilson coefficients to the structure functions F 2 (x,Q 2 ) and F L (x,Q 2 ) in analytic 
form. For F 2 (x, Q 2 ) this calculation has been performed to 2-loop order in [6,7] and for F L (x, Q 2 ) 
to 3-loop order in [8]. In the region Q 2 ^> m 2 the heavy flavor Wilson coefficients for deep- 
inelastic scattering factorize into massive operator matrix elements A^di 2 /m 2 ) and the massless 
Wilson coefficients Ck(Q 2 / ji 2 ) [9-11] for all but the power suppressed contributions. The massive 
operator matrix elements are universal and contain all the mass dependence in the logarithmic 
orders and the constant term. The process dependence is due to the massless Wilson coefficients. 
In the case of the structure function F 2 (x, Q 2 ) the asymptotic heavy flavor contributions become 
quantitatively very close to those obtained in the complete calculation [4, 12] at LO and NLO 
already for Q 2 >10 m 2 . These scales are sufficiently low and match with the region analyzed in 
deeply inelastic scattering. 

In the present paper we perform a first step towards the 3-loop heavy flavor Wilson coeffi- 
cients for the structure function F 2 (x, Q 2 ). The renormalization of the massive operator matrix 
elements to 3-loop order encounters also the contributions of 0(e) at 0(a 2 ), which have not yet 
been calculated before. 1 The 2-loop 0(e) terms form finite contributions to the 0(a^) matrix 
elements with the single pole terms emerging at 1st order. We extend the work presented previ- 
ously in Ref. [7]. For the calculation of the 0(e) 2-loop contributions our representation which 
is based on hypergeometric integrals was extended straightforwardly. However, many more in- 
finite nested sums, which contain the Mellin variable N, had to be evaluated for the first time, 
since other available techniques [15-17] could not be used for this purpose. We applied both 
suitable integral representations and the summation package Sigma [18], which solves these sums 
in IlS-fields. In the result all sums can be expressed in terms of nested harmonic sums [15,19]. 

The paper is organized as follows. In section 2 the structure of the heavy flavor contributions 
to the deeply inelastic structure function is summarized for the kinematic region Q 2 ^> m 2 . The 
renormalization of the massive operator matrix elements to 3-loop order is described in section 3. 
In section 4 the 0(e) contributions to the 2-loop operator matrix elements are calculated. Sec- 
tion 5 contains the conclusions. In the appendices we present details of the calculation, newly 
derived infinite sums and related functions depending on the Mellin parameter N, and a further 
check on our result comparing the Abelian part of the first moment with the corresponding part 

In the massless case the off-shell operator matrix elements were calculated to this order for space-like 
momenta in the MS-scheme for unpolarized and polarized deeply inelastic scattering in [13,14], which are needed 
in the calculation of the 3-loop anomalous dimensions. 
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of the on-shell photon propagator, 



2 Basic Formalism 

In the twist-2 approximation, the deep-inelastic nucleon structure functions F n (x, Q 2 ), n — 2, L, 
are described as Mellin convolutions between the parton densities fj(x,fi 2 ) and the Wilson 
coefficients C\ (x, Q 2 / fx 2 ) 

Fn(x,Q 2 ) = ^<Z(x,^®Mx,S) (1) 

to all orders in perturbation theory due to the factorization theorem. Here fi 2 denotes the 
factorization scale and the Mellin convolution is given by the integral 

[A <g> B] (x) = [ dx x I dx 2 S(x - x x x 2 ) A(x 1 )B(x 2 ) . (2) 
Jo Jo 

The distributions fj refer to massless partons and the heavy flavor effects are contained in 
the Wilson coefficients only. As was shown in Ref. [6] in the region Q 2 ^> m 2 all non-power 
contributions to the heavy quark Wilson coefficients obey 

where C^. k (Q 2 / /i 2 , x) are the Wilson coefficients for massless partons and A\ i {m 2 / ji 2 , x) are the 
massive operator matrix elements. Here \i refers to the factorization scale between the heavy 
and light contributions in O n . In the convolution (3) only those terms are accounted for which 
contribute to the respective heavy flavor Wilson coefficient functions. The index fl denotes the 
flavor-decomposition and labels the pure-singlet and gluon contributions (PS,G) and three non- 
singlet (NS 1 * 1 , NS V ) combinations. Due to the fact that the diagrams considered here contain 
one heavy quark line, at 0(a s ) only Aq g contributes. Beginning with 0(a 2 ) there is also the 
pure-singlet Aq S q and the non-singlet term Aq^ + , while at 0(a 3 s ) also the two other non-singlet 
terms Aq^ and Aq^ v contribute. The corresponding combinations of quark distributions for the 
singlet and non-singlet terms are Q 2 ) and q NS \x,Q 2 ) with 

£(*,Q 2 ) = ^[g fc (^,g 2 )-g fc (a:,Q 2 )] (4) 

k=l 

Q 2 ) = Q 2 ) ± qjx, Q 2 )} - [q n (x, Q 2 ) ± qjx, Q 2 )] (5) 

N f 

q m \x,Q 2 ) = [<Ik(x, Q 2 ) ~ Q k (x, Q 2 )] ■ (6) 



fc=i 



A^; denotes the number of light quark flavors. The massless Wilson coefficients were calculated 
in [9-11] to 3-loop orders. The massive operator matrix elements are process independent 
quantities. The factorization (3) is a consequence of the renormalization group equation. In 
Mellin space the operator matrix elements and light flavor Wilson coefficients obey the 
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following expansions 



(2 \ 00 
^■) = <w> = 5 M +E o ^!f» < = 9»» ( 7 ) 
^ ' 1=1 

^ + * = 9>9 (8) 

v/i / 1=1 

of the twist-2 flavor singlet, non-singlet and gluon operators 0^ s,s ' g between partonic states 
which are related by collinear factorization to the initial-state nucleon states \N). The local 



operators are given by 
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. . D^ N - 1 F a ^ N (z)} - Trace Terms 



(9) 
(10) 



Here S denotes the operator which symmetrizes all Lorentz-indices and = <9 W — gt a A a is the 
covariant derivative, q(z), q(z) and F a ^ v {z) denote the quark-, anti-quark field and the gluon 
field-strength operators, with g = (47ro; s ) 1 / 2 = (167r 2 a s ) 1 / 2 the strong coupling constant, t a the 
generators of SU(3) C , and A r the Gell-Mann matrices of SU(3)f- The Feynman rules for the 
operator insertions are given in [7, 20] . 



3 Renormalization of the Matrix Elements 

The massive operator matrix elements contain ultraviolet and collinear singularities which have 
to be renormalized. Charge-, mass-, operator-, and wave function renormalization have to be 
performed. Collinear singularities appear in those parts of the diagrams with vertices which link 
only to massless lines, and are specific to the particular classes of diagrams. Since in the present 
case at least one closed fermion line is massive, collinear singularities appear only at 0(a 2 ). The 
un-renormalized massive operator matrix elements read 

00 

4 = % + $>*i;f . (12) 

fc=0 

Here a s denotes the bare coupling constant. To 2-loop order, the corresponding diagrams were 
given in Ref. [6]. Here one has to distinguish one-particle irreducible and reducible diagrams, 
which both contribute in the calculation. We would like to remind the reader the background of 
this aspect. 

If one evaluates the heavy-quark Wilson coefficients in an usual Feynman- diagram calculation, 
the matrix elements are given by diagrams of the type depicted in Figure 1. The incoming gluon 
is factorized from the nucleon, i.e. we assume the parton life-time tl being much longer than the 
interaction time 77 of the virtual photon with the nucleon. As is well-known [21], this condition 
is fulfilled whenever k\ Q 2 and neither the Bjorken variable x is very small (x <jt 1) nor large 
(x 76 1). This is the case performing the Bjorken limit and applying the collinear parton model, 
in which the incoming massless partons are dealt with as on-shell particles. Also in this case, 
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Figure 1: Massive quark self-energy correction to virtual scalar-gluon scattering 



self-energy diagrams for the incoming parton lines are present. However, one may factorize these 
contributions into the non-perturbative parton densities at leading twist, resp. parton correlation 
functions at higher twist, since these contributions are virtual and are always present whatever 
hard scattering cross section is considered. They do not form a heavy quark signature which can 
be identified in a subspace of the complete final-state Fock-space emerging in deeply inelastic 
lepton-nucleon scattering. This procedure was adopted in Ref. [4]. One consequence is that at 
0(a 2 ) there are no diagrams with two fermion lines, resp. at 0(a^) none with three fermion 
lines in the general heavy flavor Wilson coefficients. The situation is different in case of the 
operator matrix elements obtained after the light-cone expansion is being performed. Here, the 
line between the two virtual photon- or weak gauge boson vertices is contracted. This line may 
contain virtual corrections, see e.g. Figure 1, which would be lost in the process of contraction. 
They have to be accounted for in attaching these self-energies to outer lines of the contracted 
diagram, see Figure 2. From the case of the fermion-fermion anomalous dimension at leading 
order these aspects are known for long [22-24]. 




Figure 2: 0(a s ) Self-energy correction due to massive quarks for the operator matrix element AgJ . 

In kinematic regions, where higher twist effects can be safely neglected [3,25] and at suffi- 
ciently high scales Q 2 , the scaling violations of deeply inelastic structure functions are due to 
the running coupling constant and heavy quark mass effects, after target mass effects [26] have 
been accounted for. We will further assume that we are in a region where power corrections 
due to heavy quarks are negligibly small, i.e., the heavy quark effects contribute logarithmically 
oc In (m 2 / jj 2 ), I > 0. In this region one may express the structure functions Fi(N, Q 2 ) in Mellin 
space by 

Fi{N, Q 2 ) = J2 C ^(N,Q 2 /^a s (fi 2 .))- [g / (iV,^/^,a s (/i r 2 ))+^(iV, / iJ/ yU 2,a s ( / i r 2 ))] 
i=i 

+C i>g (N, Q 2 /fi 2 f ; a s (/i 2 )) ■ g{N, /ij/^u 2 , a s (/i 2 )) 
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+ J2H l , q (N,Q 2 /fi 2 f ;a s (fi 2 r )) ■ [ Ql (N, fi 2 ,a s (fi 2 r )) + q l (N, )j, 2 ,a s (tf))] 

h=l 

+H ii9 (N, Q 2 /fi 2 ; a s (fi 2 r )) ■ g(N, a s {£)) . (13) 

Here fi r and /// denote the renormalization and factorization scales, respectively, no is a hadronic 
scale, and q, q and g denote the quark- and gluon distribution functions. Furthermore, the heavy 
quark Wilson coefficients factorize according to (3), which is described by the scale /i. In the 
following we identify all these scales /i — ii r — /if. Since the structure functions Fi(N,Q 2 ) do 
not depend on these scales they obey the following renormalization group equation (RGE) [27] 

[D + 2 7j ] F n (N,Q 2 )=0, (14) 

where the differential operator D is defined by 

/3(a s ) denotes the /^-function, 7 m (a s ) the mass anomalous dimension, and 7j denote the anoma- 
lous dimensions of the quark fields. Here we discuss the case of conserved currents, which have 
vanishing anomalous dimensions. 2 

The RGE for the Wilson coefficients and the parton distributions read [23, 28] 

[23^-7S(a s )]Q(7V,g 2 / M 2 ) = (16) 
[(T> + 2 lj )5 kk , + ^ k ,(a s )}f jk ,(N,v 2 /iil) = 0, (17) 

with 

f jk (N,» 2 /i4) = (j\O k \j) . (18) 
In the following we describe the renormalization to 0(a 3 s ). 



3.1 Charge Renormalization 



We perform the charge renormalization in the MS-scheme. This allows to compare the results 
obtained in the QCD analysis of deeply inelastic scattering data with analyzes of other data. 
The bare coupling constant a s is expressed by the renormalized coupling a s in the MS scheme 
by 



a s (e) = Z 2 (e,fi 2 )a s (fi 2 ) 



a s (^ 2 ) [l + 5a Syl a s (^ 2 ) + 5a Sy2 a 2 s (fi 2 )] + 0(a i s ) 



s 2 



e 2 e 



with Z g the Z-factor for the strong charge. Here the spherical factor S £ is given by 



S F 



cxp 



- { lE - 1h(4tt)) 



(19) 
(20) 

(21) 
(22) 



Calculating the evolution of the transversity structure function h\{x,Q 2 ) using the forward Compton ampli- 
tude, this is not the case, cf. [29]. 
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with 7s the Euler-Mascheroni number, e = D — 4, and D the dimension of space-time. (3q and 
Pi [30] denote the first expansion coefficients of the /3-function in the massless case 



da s {fi 2 ) 
d\n{n 2 ) 

Po 

Pi 



2 
11 



2^- 

k=0 



-Ca - ^T F n f 



34 r2 
3 A 



A^-C A + C F ^T F n } 



(23) 

(24) 
(25) 



The color factors for SU{3) C are C F = (iV c 2 - 1)/(2N C ) = 4/3, C A = N c = 3,T F = 1/2. n f 
denotes the number of active flavors. The renormalized coupling constant is obtained absorbing 
Z g into the bare coupling g. 

In [6] a slightly different point of view was taken, including mass effects in the evolution 
of a s (fj, 2 ) which usually means to choose another scheme, as e.g. the MOM-scheme [31]. To 
maintain the Slavnov- Taylor identities of QCD the calculation has to be performed using the 
background- field method [32] in [33]. 3 Since various mass scales contribute even in case power 
corrections can be disregarded, to treat a s including mass effects is also somewhat non practical. 
We treat the corresponding mass effects explicitly, which is outlined in Section 3.4 below. 



3.2 Mass Renormalization 

We choose the on-mass-shell scheme for quarks. In case of the heavy quarks the bare mass m 
is related to the renormalized mass by 



m 
Smi 

5m 2 



m + a s 5mi + a^<5m 2 + 0(a 



C F S £ m 



m 



/// 



2\ e/2 



(26) 
(27) 



C F Stm ( ^ £ 



- (18C F + 22C A - 8T F (iV, + N h )) 



1/45 91 
+ - [~Y Cf + ~2 Ca ~ UT f( N i + N h) 

iqq k\ \ 
+( ( _ - _C 2 + 481n(2)C 2 - 12 C 3 + C A 



fif)5 5 

■^ + |C 2 -241n(2)C 2 + 6C3 



45 \ /69 

v; ( Y + 10C2J +nJ y -U( 2 



(28) 



(27) is easily obtained. The pole terms to (28) were given in [35], after charge renormalization, 
and the constant term in [36], see also [37]. The 3-loop corrections were given in [38]. The 
renormalized mass is obtained absorbing Z m into the bare mass m. Heavy quark mass effects 
occur also for massless quark self-energies, see Section 3.4. 



3 Earlier calculations [26, 34] illustrated this reporting different expressions for Z g depending on the vertex 
considered. 
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3.3 Operator Renormalization 

The local operators which emerge in the light cone expansion contain ultraviolet divergences. 
These are renormalized by the following Z-factors for the flavor non-singlet (NS), singlet (S), 
and pure singlet (PS) contributions. The formulae are partly generic and have to be adapted, 
e.g. for the three flavor non-singlet contributions. Here we suppress the argument N in the 
anomalous dimensions 7^. 



Z NS (N,a s ,e) = 



Z qq (N, a s ,e) 



1 . c 7NS,0 . 2 C 2 

1 + a s S £ + a s S £ 



1 fl 



^2 ( 2 1ns >° + ^ o7NS '° ) + ^7ns,i 



2e 



+a 3 s S 3 



1 4 

^3 ( ^7ns,o + A)7ns,o + ^PolNSfi 



1 fl 



-o o7ns,o7ns,i + ~A)7ns,i + t;/3i7ns,o + tj-7ns,2 



e 2 \2 



1 + a s S £ ^ + a 2 s S 2 £ 



3e 



(29) 



(7^,0 + 7g 9 ,o7<?g,o) + A)7gg,0 



+a 3 s S 3 



4 

+ ^Phqq,0 



(llqfi + 2 lqq,0lqgfilgq,0 + 1 qgfll ggfil gqfi) + A) {llqfi + Iqgfilgqfl) 



+ 



1112 

-jlqqfilqq,! + ^Iqgfilgq,! + g7<W,l7jra,0 + g (A>7 OT ,1 + ATw.o) 



7gg,2 ' 

3e 



(30) 



Z ?9 (iV,a s ,£) = 



1 



(7g 9 ,07<y<y,0 + 7gg,o7g 9 ,o) + A)7g<?,0 



g (799,07^,0 + 7gg,07g 9 ,07<?<?,0 + 7g 9 ,07 9 g,07g 9 ,o) 

4 



+A) (7<7 9 ,07<y<y,0 + Iqqfilqgfi) + ^Polqgfi 



1 

+ ^7 



1 2 

- (7g 9 ,l7ss,0 + IqqAlqgfi + 2jqq,0jqg,l + 27g 9 , 7 9 <?,l ) + g (A>7«ff,l + ft7?S,o) 



7gg,2 • 

3e . 



(31) 



Z OT (iV, a s ,e) = a s £ £ 



7<?g,0 



+ o?5? <! ^ 



2 (7 9 g,o7gg,0 + 7<?<?,o7 9 g,o) + A>7s?,0 



+ 2£ 7 OT ,i 



« (7 OT ,07jg,0 + 7 33 ,o7 9 g,o7gg,0 + 7g< ? ,07?< 7 ,07< 79 ,o) 



+A) (7 S g,o7gg,0 + 7 ra ,o7 S g,o) + ^llgqfi 



+ - 



1 2 

g (7 S g,l7gg,0 + 7<y<y,l7 9 g,0 + 27 ra ,o7 9 g,l + 27 OT ,o7gg,l) + g (A>7ot,1 + ft7«,o) 



7gg 

' 3e 



(32) 
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Z gg (N,a s ,e) = 1 + + «»5 £ 2 2 



+ ^0 7 5 3,0 



r (7gg,0 + 2733,0739,0793,0 + 733,07<73,0733,o) + A) (iggfl + Igqfilqgfi) 



£ 2 



1112 

2733,0733,1 + 3739,07^3,1 + g739,l7?3,0 + 3 [Polgg,! + Pllggfl) 



7gg,2 • 
3e 



(33) 



The pure-singlet operator has the following Z-f actor. 

1 1 



PS 

-^Iqgfilgqfl + 7^7 g9 ,l 



799,2 



(2799,0793,0739,0 + 793,0733,0739,o) + (2793,0739,1 + 793,l7 9 9,o) 



3e 



(34) 



The anomalous dimensions 7ij,fc(iV) are related to the splitting functions by 

7«,*W = - fdzz N -'P^{z). (35) 
Jo 



The renormalized operators are obtained absorbing Z^s, Zj^g, an d 7 m to the bare operators, 
resp. operator matrix elements. 



3.4 Wave Function Renormalization 

The external legs of the operator matrix elements are treated on-shell to be able to apply their 
factorization from the nucleon wave-functions in the light cone expansion as outlined above. 
Here the mass scale is set by a heavy quark mass. To the operator matrix elements also one- 
particle reducible diagrams contribute. If either the self-energy insertion on the external legs or 
the remainder diagram contain only massless lines, with the exception of the tree-level terms, 
the diagrams are vanishing since one of the factors has no scale. I.e. finite contributions are due 
to the self-energy insertions containing a massive line. The corresponding corrections are due to 
the massive contributions to the massless quark self-energy up to 3-loop order and the gluon 
self-energies up to 2-loop order. The former terms emerge in case of the flavor non-singlet terms 
A^qQ® and the latter in Aq 1 , while the pure singlet contributions A^'^ obtain no corrections. 



3.4.1 Massless External Quark Lines 

The 2-loop correction reads 



Sg } = a%,T F C F [ — ) S'i { - + ^ + 



2 \ £ 

f< 2 ' 



6 



89 
72 



2 



e + 0(e 2 ) \-ijf. 



(36) 



At 0{a 2 s ) this contribution implies, that the 1st moment of the non-singlet operator matrix 
element vanishes. 
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3.4.2 External Gluon Lines 



The gluon vacuum polarization is given by 

n£(<7) = [-g,uq 2 + q,qu}H ab (q 2 ) , 
with a and b the color indices. The 1-loop and 2-loop corrections read 



nft(o) 



(37) 



(38) 



ng)(o) 



-ia 2 s 5 ab T F S l 



2\ £ 



12 13 /35 

T + y + i2 +3C2,: 



e 2 e 



13 
12 



-C2 " 



169 5 Cs , 
144 + 4 C2 _ 3 ' 



(39) 



The C^-term can be compared with a corresponding contribution in the photon propagator, 
(180), before mass renormalization. At 2-loop order the diagrams u and v from [6] and the term 
Zqg'^A^oo combine to 



*99 



m 



with 



Po,h — ' 



(40) 



(41) 



(40) yields a finite contribution oc Tj. in the MS scheme. Our treatment differs from that in 
Ref. [6] as we do not include the mass effects of (40) into the running coupling, because we have 
chosen to define it in the MS-scheme. This is convenient for direct comparisons of the parton 
densities and the QCD-scale Aq CD measured in other analyzes of hard scattering cross sections. 



with 
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3.5 Mass Factorization 

The mass singularities are factored into the functions T NS , T^ s an d r g <?,ps, respectively. If all 
quarks were massless these functions were given by 

(42) 
(43) 
(44) 

(45) 
(46) 
(47) 



r 9g ,ps = z gg ] PS , 

7ns,o 



7-1 



T m (N,a s ,e) = 1 - a s S £ ^ + a 2 s S 2 

e 



Fij,s(N, a s ,e) 
F qq ,ps{N, a s ,e) 



Oij ds^e - r o s o £ 



' 1 (\ 



e 2 \2 



-a 2 s S 2 



1 1 

-^Iqgfllgqfl + 7^7<7<?,PS,1 
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In the present calculation at least one quark line is massive in each diagram. Therefore the 
T-matrices (45-47) apply to the parts of the diagrams which contain massless lines only, which 
are at most 2-loop sub-graphs. The mass factorization is therefore different in various sub- 
classes of contributing Feynman diagrams. The functions r^s, T^g, an d T^ps do thus enter 
the renormalization of the operator matrix elements only in products with other functions. The 
singularities contained in Tns, r^g, and r^pg are absorbed into the bare parton densities, which 
become scale-dependent in this way. 



3.6 The renormalized operator matrix elements 

The operator matrix element reads after charge and mass renormalization 



5ij + a s A\f + a 2 s 



'dm 



2 dm 2 



(48) 



The renormalized operator matrix elements are obtained removing the ultraviolet singularities 
and collinear singularities of the operator matrix elements, 



Aj = Z^A k ^=5 lJ +a s A < ff + a^A < ff+^ 



(49) 



Here self energy insertions containing massive lines in the external legs of the operator matrix 
elements have to be kept. 



4 The O(e) Contributions 

The 0(e) contributions to Aq g , Ag S q and A^ Q at 0(a 2 s ) contribute to these quantities at 0(a 3 s ) 
in combination with the various single pole terms emerging at 1-loop, as outlined in Section 3. 
The diagrams to be evaluated are shown in [6,7]. The results for the individual un- renormalized 
diagrams in 0(e) are given in Appendix A. As outlined before in Ref. [7], we calculate the 
massive operator matrix elements performing the Feynman-parameter integrals directly, i.e., 
without using the integration-by-parts method [39] which was applied in [6] up to the terms 0(e°) 
before. We obtain representations in terms of generalized hypergeometric functions [40], which 
may be expanded to the desired order in e. With increasing depth in e, more and more involved 
nested infinite sums are obtained, which depend on the Mellin-parameter N. These sums can be 
summed applying analytic methods, as integral representations, and general summation methods, 
as encoded in the Sigma package [18]. We applied both methods to evaluate the sums which 
emerge at 0(e). The underlying algorithms of Sigma are based on a refined version [41] of Karr's 
difference field theory of ITE-fields [42]. In this algebraic setting one can represent completely 
algorithmically indefinite nested sums and products without introducing any algebraic relations 
between them. Note that this general class of sum expressions covers as special cases, e.g., the 
harmonic sums [15,19] or generalized nested harmonic sums cf. [43-46]. Given such an optimal 
representation, by introducing as less sums as possible, various summation principles are available 
in Sigma. In this article we applied the following strategy which has been generalized from the 
the hypergeometric case [47] to the ITE-field setting. 
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1. Given a definite sum that involves an extra parameter N, for typical sums see the Ap- 
pendix B. We compute a recurrence relation in N that is fulfilled by the input sum. The 
underlying difference field algorithms exploit Zeilberger's creative telescoping principle [47] . 

2. Then we solve the derived recurrence in terms of the so-called d'Alembertian solutions [47]. 
Since this class covers the harmonic sums, we find all solutions in terms of harmonic sums. 

3. Taking the initial values of the original input sum, we can combine the solutions found 
from step 2 in order to arrive at a closed form representation in terms of harmonic sums. 

A detailed example for the sum (120) with all its computation steps has been carried out in [48]. 
In Appendix B we present the details for the calculation of a further example, Eq. (151), in 
section (B.8). 

The results for new sums contributing are listed in Appendix B. In the calculation also more 
well-known sums are occurring which were found before in [7] or can be easily solved using the 
FORM-code [49] summer [15]. 

The 0(e) contribution to Aq q reads : 

agW = T F cJ N m + wx1 2) (i6g 2 ,i,i - as 8 ,i - s^fr + 3S 4 - is 3 s 1 - \sl - S 2 S\ 




-S, + 2( 2 S 2 - 2CA - -C3S1J - 8 N 2 {N+1){N + 2) S ^ + 3Ni(N + 2) S 



2 37V 4 + 48iV 3 + 43iV 2 - 22iV - 8 3iV + 2 Sj_ 

+ 3 N 2 {N + 1) 2 {N + 2) 3+ N 2 (N + 2) 2 1+ N 2 ^ 2 

2 (N 2 + N + 2) (3N 2 + 3N + 2) I\ 

+ 3 N 2 (N + 1) 2 (N + 2) Cs + N 3 (N + l) 3 (iV + 2) 

N 4 - 5iV 3 - 32iV 2 - 18iV - 4 2 5iV 6 + 15iV 5 + 36iV 4 + 51iV 3 + 25iV 2 + 8iV + 4 
+ N 2 (N + 1) 2 (N + 2) 1 iV 3 (iV + l) 3 (iV + 2) ^ 2 

2N 5 - 2N 4 - UN 3 - 19iV 2 - 44iV - 12 P 2 



+T F C A 



N 2 (N + 1) 3 (N + 2) N 5 (N + 1) 5 (N + 2) 

N 2 + N + 2 



N(N + l)(N + 2) 

2 



x \16S. 2 , 1A - 4S- 2)lil - 85_ 3) i - 8,S_ 2 , 2 - 4S 3 ,i - + 9S 4 
-16^2,!^ + ^S 1 S 3 + 40" Si - 8/3% + ^Sl - 8P'S 2 + 5S 2 S 2 + ^ - ^-S 1 ( 3 

-2s 2(2 - 2s& - - ^ci) + { T + \n N ^) 2 + ? - 4 ^) 

16 N 5 + ION 4 + 9N 3 + 3iV 2 + 7N + 6 3iV 3 - 12iV 2 - 27iV -2 

-S 3 + 2 — Ar/Ar ; lN9 / Ar , ^ 2 — ^2^1 



(N — 1)N 2 (N + 1) 2 (N + 2) 2 3 N(N + l) 2 (iV + 2) 



2iV 3 + 8iV 2 + lliV + 2^ „ N 2 + N 



3 N(N + l) 2 (iV + 2) 2 1_ (N+ l) 2 (iV + 2) 2 ^ 2 1 

2 9N 5 - ION 4 - UN 3 + 68iV 2 + 24iV + 16 N 4 + 2iV 3 + 7iV 2 + 22iV + 20 , 

3 (N — 1)N 2 (N + 1) 2 (N + 2) 2 ^ 3 + 8 (N + l) 3 (iV + 2) 3 " 
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P 3 2P 4 



(iV- l)iV 3 (iV+ l) 3 (iV + 2) 3 z (iV- l)iV 3 (iV + l) 3 (iV + 2) 



7V(7V + l) 3 (iV + 2) 3 1 N(N +1) 4 (N + 2) 
2P 7 



(N — 1)N 5 (N + 1) 5 (N + 2) 5 
Here the argument TV of the harmonic sums, and (N + 1) in the function 



(50) 



(51) 

S-i(N) = (-l) N p(N + l)-\n(2) (52) 
and in the polynomials Pi(N) was omitted as well as the factor 



, /m 2 \ e 



In the gluon and pure-singlet case we did not write the overall factor 

i + (-ir 

2 

It does not emerge generically in the non-singlet case. In accordance with the light-cone ex- 
pansion only even integer moments contribute in the present case. The polynomials in Eq. (50) 
are 

P 1 = 3iV 6 + 30iV 5 + 15iV 4 - 64iV 3 - 56iV 2 - 20N - 8 , (54) 
P 2 = 2AN 10 + 136iV 9 + 395iV 8 + 704iV 7 + 739iV 6 + 407iV 5 + 87iV 4 + 27iV 3 + 45iV 2 + 24iV + 4 , 
P 3 = AT 9 + 21iV 8 + 85iV 7 + 105iV 6 + 42iV 5 + 290iV 4 + 600iV 3 + 456iV 2 + 256iV + 64 (55) 
P 4 = (jV 3 + 3iV 2 + 12iV + 4)(iV 5 -iV 4 + 5iV 2 + iV + 2) , (56) 
P 5 = N 6 + 6N 5 + 7N 4 + AN 3 + 18N 2 + 16 N - 8 , (57) 
P 6 = 2N 8 + 22N 7 + 117N 6 + 386N 5 + 759N A + 810N 3 + 39QN 2 + 72N + 32 , (58) 
P 7 = AN 15 + 50N U + 2Q7N 13 + 7Q5N 12 + 1183N 11 + Q82N 10 -82QN 9 - 1858N 8 (59) 
-1116iV 7 + 457iV 6 + 1500iV 5 + 2268iV 4 + 2400iV 3 + 1392iV 2 + 448iV + 64 . (60) 

The flavor non-singlet and pure-singlet contributions read : 

^8,(2) _ Tr / 4 9 + 4 9 ^ 8 SC 2 ° 9 20 3^+3^+2 112 

3N 4 + QN 3 + A7N 2 + 20N -12 _ 656_ c? P 8 \ 

+ 18iV 2 (iV + l) 2 C2 ~ ~81 1 + 6A8N*{N + l) 4 j ' ( ' 

P 8 = 1551iV 8 + 6204iV 7 + 15338iV 6 + 17868iV 5 + 8319iV 4 

+944iV 3 + 528iV 2 - 144iV - 432 . (62) 
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T F C F l -2 



(5iV 3 + 7N 2 + 4N + 4)(iV 2 + 5N + 2) 



(iV-l)iV 3 (iV + l) 3 (iV + 2) 2 
(iV 2 + iV + 2) 2 



(2S 2 + C2) 



3 (AT — l)iV 2 (iV + l) 2 (iV + 2) 



( 3 * 3 + C3 ) +2 7iv3I 



(iV-l)iV 5 (iV+ l) 5 (iV + 2) 4 



P 9 = 5iV n + 62iV 10 + 252iV 9 + 374iV 8 - 400iV 6 + 38iV 7 - 473iV E 
-682iV 4 - 904iV 3 - 592iV 2 - 208iV - 32 . 



(63) 
(64) 



The harmonic sums contributing to the individual diagrams, see Appendix A, are listed in 
Table 1. 

Table 1: Complexity of the results in Mellin space for the individual diagrams in the unpolarized case, 
cf. [7], up to 0(e) 



Diagram 


£1 


£2 


£3 


C 
04 


£-2 


£-3 


C 

Q—4 


£2,1 


^-2,1 


£-2,2 


£3,1 


£-3,1 


£2,1,1 


£-2,1,1 


A 




+ 


+ 
























B 


+ 


+ 


+ 


+ 








+ 






+ 




+ 




C 




+ 


+ 
























1 J 


1 


1 
\ 


1 
\ 










1 
\ 














E 


+ 


+ 


+ 










+ 














F 


+ 


+ 


+ 


+ 








+ 










+ 




G 


+ 


+ 


+ 










+ 














H 


+ 


+ 


+ 










+ 














I 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


J 




+ 


+ 
























K 




+ 


+ 
























L 


+ 


+ 


+ 


+ 








+ 






+ 




+ 




M 




+ 


+ 
























N 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 





+ 


+ 


+ 


+ 








+ 






+ 




+ 




P 


+ 


+ 


+ 


+ 








+ 






+ 




+ 




s 




+ 


+ 
























T 




+ 


+ 
























PS a 




+ 


+ 
























PS 6 




+ 


+ 
























NS a 






























NS fe 


+ 


+ 


+ 


+ 























Here we have already made use of the algebraic relations [50]. Moreover, two of the sums, 
S-2,2(N) and S^±(N), can be related by structural relations [51] to other harmonic sums, i.e., 
they lie in corresponding equivalence classes and may be obtained by either rational argument 
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relations and/or differentiation w.r.t. N. Reference to these equivalence classes is useful since 
the representation of these sums for iV e C needs not to be derived newly, except of differentiation 
which is easily carried out. Therefore the two-loop massive operator matrix elements to 0(e) 
depend on six basic harmonic sums. 



diagram/order 


l/e 2 


1/e 


1 


£ 


e 2 


B 


N = 2 


o 

— 


4.00000 


—8.82690 


2.4 1/28 


— o.Dyo2o 




N = 6 


-7.73333 


0.81936 


-8.89777 


-1.84111 


-7.25674 


C 


N = 2 


—4 


la. 8 


9 fil 71 C 

—3.61 1 15 


17 Q 9 1 n Q 

1 1 .33108 


O.ZD440 




N = 6 


-1.33333 


8.26984 


-1.34024 


7.12612 


1.38782 


D 


N = 2 


o 

— 5 


1 .80006 


C O A C A O 

—6.34542 


A 71 

4. 1 123b 


— 2.18586 




N = 6 


-2.66666 


-0.69523 


-2.60657 


-1.74990 


-2.37611 


E 


N = 2 


5.5588a 


— 11.2593 


9.82824 


— 12.8921 


O 9A1 A C 

2.39145 




N = 6 


2.93878 


-4.24257 


3.39094 


-4.3892 


0.826978 


F 


N = 2 


2.66667 


-4.88889 


9.1707 


-1.2618 


8.10105 




N = 6 


1.65714 


-3.43645 


6.12836 


-0.81895 


5.43478 


G 


N = 2 


2.66666 


-9.55555 


4.59662 


-8.92015 


1.07313 




N = 6 


0.57142 


-2.00204 


1.04814 


-1.89142 


0.32219 


II 


N = 2 


-2.66666 


3.55555 


-3.76328 


4.57782 


-2.02934 




N = 6 


-0.57143 


0.83061 


-0.9478 


1.08896 


-0.55443 


I 


N = 2 



















N = 6 


0.34285 C A 
+0.C F 


-0.38979 C' A 
+0.30793 C F 


0.45862 C A 
-0.00949 C F 


-0.33621 C' A 
+0.32243 C F 


0.27344 C A 
+0.10405 C F 


L 


N = 2 


22 22222 


—8.07407 


21 63853 


—1 09400 


14 38017 

It:. 'J U y > J. 1 




N = 6 


9.42675 


0.69179 


9.99619 


4.12879 


8.49503 


M 


N = 2 


-0.44444 


-1.66666 


-0.25375 


-2.19775 


-0.72592 




N = 6 


0.09977 


-0.21941 


0.12862 


-0.23387 


0.03250 


N 


N = 2 


5.77778 


-1.70370 


4.24641 


-0.06679 


2.27892 




N = 6 


4.26462 


0.17238 


4.61023 


1.71768 


3.87767 





N = 2 


-4.44444 


3.14815 


-4.96659 


2.01274 


-3.0498 




N = 6 


-1.03605 


0.56098 


-1.12332 


0.27406 


-0.73034 



Table 2: Numerical values for the moments N = 2,Q for the the expansion of the un-renormalized 
matrix element A 2 Qg for the terms 0(1/ e 2 ) to 0(e 2 ) for individual diagrams. 

In the complex Mellin-iV plane these functions, up to more simple terms due to the soft- and 
virtual corrections, are meromorphic functions with poles at the non-positive integers, which 
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possess both an analytic regular asymptotic representation and recursion relations, through 
which they may be calculated. To apply the results obtained in the present calculation in Mellin 
space in a QCD-analysis of deeply inelastic structure functions, their scale evolution is first 
evaluated in Mellin space, incorporating the heavy flavor Wilson coefficients for complex values 
of N. This requires analytic continuations of the corresponding harmonic sums as worked out 
in [52]. The result in x-space is obtained by a single, fast numeric Mellin-inversion performed 
by a contour integral around the singularities of the problem located at the real axis left to some 
value r. In case one wants to include small-x resummations as well, this can be done in a similar 
way, see [53]. 

We performed an independent check on our calculation evaluating fixed moments in N for the 
un-renormalized diagrams using the Mellin-Barnes method [54,55]. Here we use an extension 
of a method developed for massless propagators in [56] to massive on-shell operator matrix 
elements [57]. The Mellin-Barnes integrals are evaluated numerically using the package MB [58]. 
In Table 2, we present the moments N = 2 and N = 6 for the more difficult two-loop diagrams, 
cf. [7], for the 0(1 /e 2 ) to the 0(e 2 ) terms. 

A further test for the Abelian part of the first moment of the un-renormalized massive oper- 
ator matrix element Aq q after mass renormalization, i.e., the term oc TfCf, can be performed 
after analytic continuation of even values of N. This term is related to a corresponding contri- 
bution of the on-shell photon polarization function as noted in [6] . We apply this method to the 
0(e) term in Appendix C and find agreement. 

5 Conclusions 

We calculated the 0(e) contributions to the massive operator matrix elements at 0(a 2 s ) which 
contribute to the heavy flavor Wilson coefficients in deeply inelastic scattering to the non power- 
suppressed contributions. In the renormalization of the heavy flavor Wilson coefficients to 3- 
loop order they contribute together with with the single pole terms at 0(a s ). These terms, and 
the 0(a 2 e) contributions to the operator matrix element A gg (N) to be published soon, form 
all but the constant terms of the 3-loop heavy flavor unpolarized operator matrix elements 
needed to describe the 3-loop heavy flavor Wilson coefficients, together with the known 3-loop 
massless Wilson coefficients [11], in the region Q 2 ^> m 2 . In the calculation, we made use of 
the representation of the Feynman-parameter integrals in terms of generalized hypergeometric 
functions in a direct calculation, without applying the integration-by-parts method. The e- 
expansion leads to new infinite sums which had to be solved by analytic and advanced algebraic 
methods. We checked our results for finite values of N, using the Mellin-Barnes method, and 
for a series of diagrams by a second program. Here, the calculation can be extended to higher 
corrections in e. For N — 1, one may compare in addition the terms oc TfCf in o!q g with 
the corresponding contribution in the 2-loop on-shell photon propagator. The terms aP\N) 
can be expressed in terms of polynomials of the basic nested harmonic sums up to weight 
w = 4 and derivatives thereof. They belong to the complexity-class of the general two-loop 
Wilson coefficients or hard scattering cross sections in massless QED and QCD found for space- 
and time-like unpolarized and polarized anomalous dimensions, massless Wilson coefficients for 
deeply inelastic scattering, parton fragmentation, the Drell-Yan process, hadronic Higgs- and 
pseudoscalar Higgs production in the heavy mass limit as well as the soft- and virtual contribution 
to Bhabha-scattering, cf. [60], and are described by six basic functions and their derivatives in 
Mellin space. Their analytic continuation to complex values of N is known in explicit form. 



16 



The package Sigma [18] proved to be a useful tool to solve the sums occurring in the present 
problem and was extended accordingly. One advantage to seek for solutions of the recurrences 
emerging in ITE-fields consists in finding irreducible structures for the representation. In the 
present calculation these were nested harmonic sums. In even more complicated single scale 
problems in higher orders, this needs not to be the case. The new basis elements, however, 
would be uniquely found applying the present procedure. 
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A The O(e) Terms for the Individual Diagrams 



In the following we list the results for the individual diagrams for comparisons and to illustrate 
the analytic structures emerging in the calculation. The calculation is performed in Feynman 
gauge. Again we suppress the factor (53) and the argument TV in the sums and polynomials. 

2iV 3 - TV - 2 /„ „ u \ 2A 



-jQg 



p 2 
P3 

A? 9 



P 5 



3TV 2 (TV + 1) TV 3 (TV+ l) 2 (TV + 2) 



2S 2 + ( 2 



TV 5 (TV + 1) 4 (TV + 2) 3 



Pi = 



-jQg 



TpC p ' 

2TV 9 - 16iV 8 - 89iV 7 - 166iV 6 - 135iV 5 - 6TV 4 + 85TV 3 + 94iV 2 + 44TV + 8 . 
TfCf< 



,(65) 



"| TV O^ 2 ' 1 ' 1 ~~ ^3,1 ~~ 8<S , 2j i<S'i — -S 3 Si + -Sf — S 2 S 2 — -Si + 6S 2 (2 



-2Slt 2 - l^xCs + \Cz) + 



TV 2 + 7TV + 2 
TV(TV + l)(TV + 2) 



(sS 2 , 1 + 2S 2 S 1 + ^S 3 1 +AS 1 C 2 ) 



TV 5 + 3TV 4 + 19TV 3 + 37TV 2 + 167V + 4 
TV 2 (TV + l) 2 (TV + 2) 2 
TV 2 + 5TV + 2 4 TV 2 

" 8 iVXiVTT)(iVT2) + 

I6P3 



So -A 



TV + 9TV 2 + 8TV + 4 



17iV + 2 



-S 3 + 



N 2 {N + 2) 2 
16P 2 



3TV(TV + l)(TV + 2) * N 2 (N + 1) 3 (TV + 2) 



:S\ 



TV(TV + l) 4 (TV + 2) 3 
TV 7 + UN 6 + 65iV 5 + 153iV 4 + 197TV 3 



134iV 2 + 44TV + 8 , 



2N 7 + 27TV 6 + 130TV 5 + 306TV 4 + 385TV 3 + 266TV 2 + lOOiV + 16 . 

7iV 3 + 29TV 2 + 15iV + 2 



T F C F ^(5S 3 -±( 3 )-2- 



N 2 {N + l)(N + 2) 



-So 



Pa = 



-jQg 



13TV 4 + 82TV 3 + 82TV 2 + TV - 6 P 4 
+ TV 2 (TV+ l)(TV + 2)(TV + 3) ^ 2 + TV 4 (TV+ l) 3 (TV + 2) 3 (TV + 3) 

32N 10 + 448TV 9 + 2177iV 8 + 5123TV 7 + 6312TV 6 + 3863TV 5 + 902TV 4 
+9iV 3 - 74TV 2 - 68iV - 24 . 

fl / 2 1 4 

— y-ASo.,! - -S3 - S 2 Si - -S l - 2Si( 2 - -( 3 



+ 



iV 4 + 8TV 3 + 43TV 2 + 36TV + 12 
TV 2 (TV + l) 2 (iV + 2) 

4P 5 



(S 2 + Sf) + 2 



TV 3 + ION 2 + 59TV + 42 
7V(7V + 1)(N + 2)(N + 3) 



C2 



4P fi 



7V 2 (7V + l) 3 (7V + 2) 2 1 7V(7V + l) 4 (7V + 2) 3 (7V + 3) 
TV 6 + 87V 5 + 797V 4 + 20 77V 3 + 2057V 2 + 967V + 20 , 

27V 8 + 247V 7 + 2627V 6 + 13717V 5 + 35147V 4 + 47757V 3 + 35447V 2 + 14047V + 240 



(66) 



(67) 



(68) 



C F - 



-2 



C, 1 „ 7V + 2 

7V(7V + 1) 

1 TV 3 - TV 2 - 87V - 36 
'3 7V(7V + l)(7V + 2)(7V + 3) 



2S 2j i + Si( 2 ) — 



5, 



3S 2 Si + S 3 ) + 



2 137V 4 + 607V 3 + 1117V 2 + 47V - 36 

3 7V(7V+ l) 2 (7V + 2)(7V + 3) 

4 TV + 3 37V 3 + 97V 2 + 127V + 4 

S3 - 2 — Ar/Ar , lW „ ; x — C2 



3(7V + 1) ; 



7V(7V + l) 3 (7V + 2) 
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9 



A 4N 5 + 1 IN 4 + 15N 3 - 86N 2 - 92N - 24 2 

+ N 2 (N + l) 3 (iV + 2)(iV + 3) 2 + N 2 (N + l) 2 (iV + 2)(N + 3) 1 



iV 2 (iV + l) 3 (iV + 2) 2 (iV + 3) 1 iV 3 (iV + l) 5 (iV + 2) 3 (iV + 3) 

P 7 = 20iV 6 + 119A r5 + 290iV 4 + 105iV 3 - 290iV 2 -212iV-24 , 
A = 8Af 7 + 62N 6 + 181iV 5 + 127iV 4 - 226iV 3 -404iV 2 - 296iV-96 , 
A = 38iV 10 + 394iV 9 + 1775iV 8 + 4358iV 7 + 6323iV 6 + 5788iV 5 + 3626iV 4 + 1462iV 3 
+100iV 2 - 184iV - 48 . 



A Q S 9 = Tj 



n C A 

c F - — 



■ (282,1,1 — 25*4 — S 2 — C2»S'2^ + 



N\ * ^ v 7 3(iV+ l)(7V + 2) 



2iV 2 - 5iV - 2 4 26iV 2 - 33iV - 10 

N 2 (N + l)(N + 2) 2,1 + 3iV 2 (iV+l)(iV + 2) 3 
2 2iV 2 -3iV + 2 / x 5AT + 2 

+ 3 iV 2 (iV + 1)(7V + 2) l^ 2 ^ 1 + N + V(iV + 1)(JV + 2) b 
2N 3 + 15iV 2 + 12N - 4 2 14iV 3 + 85iV 2 + 132iV + 52 
iV(iV + l) 2 (iV + 2) 2 1 + iV(iV + l) 2 (iV + 2) 2 ' 
2N 5 + 18N 4 + 25iV 3 + 10iV 2 + 44iV + 56 8( 2 

+ 4 AT/AT , i\q/A7- , 1 ~~ 



N(N+ l) 3 (iV + 2) 3 1 (iV + l) 2 (iV + 2) 

6iV 5 + 52iV 4 + 178iV 3 + 309iV 2 + 264iV + 84 1 
~ 8 (N + l) 4 (iV + 2) 3 J ' ^ ' 

-rQg „ n J 1252,1 - 585 3 + 3S 2 Si + Sf + 6C 2 5x + 8(3 , (9N + 8)(5iV 2 + 9N - 1) 



3(iV+ l)(iV + 2) N{N + 1) 2 {N + 2f 



37V 3 + 31iV 2 + 45iV + 8 2 17/V 2 + 47iV + 28 



iV(/V + l) 2 (iV + 2) 2 1 (iV+ l) 2 (iV + 2) 

6iV 5 + 104/V 4 + 376iV 3 + 514/V 2 + 277/V + 48 2 Ao 

+ /V(yV + l) 3 (/V + 2) 3 1 ~ (N + 1) 4 (/V + 2) 4 ' ' ! ' ' ! 

Ao = 74iV 6 + 722iV 5 + 2697/V 4 + 4960/V 3 + 4700/V 2 + 2143/V + 368 . 
C A 



4 0s - T 

h — J J 



_ f 4(iV + 3) / x 4 7V 2 -267V + 9 

2j|/V(/V + l)(/V + 2)\ 2,1 3/V(/V + l)(/V + 2)(/V + 3)' 



3iV(iV + l)(/V + 2)(iV + 3) V V 3(/V+ l)(iV + 2) 

iV 2 + 7iV + 8 iV 4 + 115iV 3 + 3357V 2 + 265iV + 84 

~ (iV + l) 2 (iV + 2) 2C2 iV(iV + l) 2 (iV + 2) 2 (iV + 3) 2 

7iV 4 + iV 3 - 31iV 2 - 37iV - 36 2 
+ N(N + l) 2 (iV + 2) 2 (iV + 3) 1 

2Al ft + > , (72) 



iV(iV + l) 3 (iV + 2) 3 (iV + 3) (N + l) 4 (iV + 2) 4 (iV + 3) 
Ai = 10iV 6 + 38iV 5 - 43iV 4 - 533iV 3 - 1267iV 2 - 1269iV - 456 , 
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367V 7 + 3447V 6 + 128 77V 5 + 21437V 4 + 818iV 3 - 2153iV 2 - 2795iV - 960 . 

2 \ TV 2 — TV — 4 / \ 

+2S_ 2 S 2 + 2S_ 2 S 1 + S^ 2 ( 2 j + 8 ^ — -j^jj — ^ 2 [2S-2,i - S- 3 - 2S^ 2 S 1 j 

1 f 327V 4- 27 
V + l)(iV + 2) ( 4(4iV + 7)Siu " 4{eN + 7)S " " 4(4Ar + 5)S ' ilSl + 2 

+| (36j v + 35) S 3S, - l -^Sl + l^+I^S? - ±5} 



3 * v y ^ i iV(JV + l) 2 (iV + 2) 



2 



27V 4 + TV 3 — TV 2 + 97V + 8 2 18iV 4 + 3iV 3 - 677V 2 - 397V + 8 

+4 T77T7 , TToTTT , ni o ^2,1 — 7T at/at , TToTTT , „\ 2 "I 



7V(7V + l) 2 (N + 2) 2 ' 3 7V(7V + 1) 2 (7V + 2) 
87V 4 - 37V 3 - 477V 2 - 29N + 8 1 37V 3 + 7N 2 - 37V - 8 3 
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7V(7V + l) 2 (7V + 2) 2 3 7V(7V + l) 2 (7V + 2) 2 

4(7V + 4)C 3 , o 7V 4 + 27V 3 + 77V 2 + 227V + 20 c P 13 S 2 

+ 8 — t D_2 + 



3(7V+ l)(7V + 2) 2 (7V+ l) 2 (7V + 2) 3 ~ z 7V(7V + 1) 3 (7V + 2) 3 

67V 5 + 407V 4 + 927V 3 + 947V 2 + 497V + 16 ~ 27V 3 + 87V 2 + 197V + 16 



7V(7V + l) 3 (7V + 2) 3 1 (7V + l) 2 (7V + 2) 

2-Pi4 c 2P 15 

+ at/at , 7TT7T7 , 7777 ^1 - 



2 



7V(7V + 1) 4 (7V + 2) 4 (7V+ l) 4 (7V + 2) 5 
+T F C F | (iv + i) 1 (iv + 2) (-32ff 2 , M + 1653,1 - 675 4 + + WS^S, 

+S 2 2 + 2S 2 Sl + \st - AS 2 ( 2 + 4^ 2 C 2 + 32C 2 + 22 4 ) - le ^^f 

4 3iV - 2 (o SS , 16 37V+1 

3 7V(7V+ l)(7V + 2) V 2 1 V 3 7V(7V + 1)(7V + 2) 3 
27V 2 + 27V-1 37V 2 + 37V + 1 

7V(7V + l) 2 (7V + 2) 2 7V(7V + l) 2 (7V + 2) 1 

1 17V 3 + 337V 2 + 387V + 14 1 
8 7V(7V + l) 3 (7V + 2) 61 J ' {<6) 

47V 6 + 107V 5 + 127V 4 + 87V 3 - 187V 2 - 377V - 16 , 

107V 7 + 1147V 6 + 5337V 5 + 13747V 4 + 21447V 3 + 20277V 2 + 10577V + 224 , 
207V 7 + 2367V 6 + 12027V 5 + 33847V 4 + 56887V 3 + 57207V 2 + 31957V + 768 . 

TC f 2 47V 2 + 47V-5 / \ Ag 

F A \ 3 7V 2 (7V + 1) 2 V 3 7V 5 (7V+ l) 5 (7V + 2) 3 

47V 5 + 227V 4 + 117V 3 + 137V 2 + 357V + 10 / \1 
+ NHN + 1)HN + 2 ) ( 2% + «■)[' < 74 > 
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28TV 10 + 148iV 9 + 342TV 8 + 285iV 7 - 212iV 6 - 114iV 5 + 1117iV 4 + 1587iV 3 + 826iV 2 
+260TV + 40 . 

TC [2 3TV 2 -23TV-20 / x At 

F A \ 3 (TV - 1)N(N + 1) 2 (TV + 2) V 3 ^ (TV - 1)N(N + 1) 5 (TV + 2) 4 

10TV 4 + 7TV 3 + 51TV 2 + 172TV + 112 



■ fa + C2) J , (75) 



(TV- l)N(N+ l) 3 (TV + 2) 2 
14TV 8 + 70TV 7 + 96TV 6 - 375iV 5 - 1493iV 4 - 1056TV 3 + 2392iV 2 + 4192TV + 1792 . (76) 

TfCa { iv ( -4S ' 2 ' i > 1 + 453,1 + \ Si + AS2 ' iSi + l SsSi + I s " + \ S2S ^ + 82(2 + 

2 6iV 3 + 5TV 2 - 4TV - 6 8iV 5 + 18iV 4 + lliV 3 + TV 2 + 6iV + 4 

~3 TV 2 (TV + l) 2 3 + iV 3 (iV+ l) 3 2 

(TV + 2)(2TV+1) 2 7TV 3 + 15TV 2 + 7TV + 4 2 2TV 3 + 5iV 2 + AN + 2 

-^1 + 2 9/ ,» 01 + 9 , , 9 C3 



TV 2 (TV + 1) 2 1 TV 2 (TV + 1) 3 3 TV 2 (TV+1) 2 

-P18C2 -P19 



N 3 (N + 1) 3 (TV + 2) TV 5 (TV + l) 5 (TV + 2) 

AN 6 + 30TV 5 + 55TV 4 + 38iV 3 + 4iV 2 - 107V - 4 , 

16iV 10 + 152iV 9 + 454TV 8 + 628TV 7 + 447TV 6 + 180TV 5 + 52TV 4 

-10TV 3 - 28iV 2 - 187V - 4 , 

_ /2TV 2 -2TV-2/ \ 27V 5 + 117V 4 + 12iV 3 + 2N 2 + 6TV + 4 / \ 

TfCa \3 N 2 (N+1) 2 [ 3S * + ty 7V 3 (7V + l) 3 (7V + 2) I 2 + C V 

67V 8 + 287V 7 + 537V 6 + 307V 5 - 147V 4 + 27V 3 + 187V 2 + 147V + 4 1 

7V 5 (7V+ l) 5 (7V + 2) J 



(77) 



(78) 



T f Ca 



^ TV (TV + 1) (~ ^~ 2 > 1 > 1 + 2S'_2,2 _ 5'_4 + 4S'_2,i5'i — 2S , _ 3 S' 1 — 2S-2S2 



-2S_ 2 Sl - S^) + 2 ^ 2 + 3 ^ + 2 : f St + 12S 2 C 2 + I6C3S1) 
1 S V 24TV(TV + l)(TV + 2) V 1 lSi S3 V 

6TV 2 + 7TV-6 / _ \ 8TV 2 + 9TV- 10 54TV 2 + 97TV -10_ 

TV(TV + 1)(TV + 2) V 2,1,1 ~ 2,1 V + TV(TV + l)(TV + 2) 3,1 4TV(TV + 1)(TV + 2) 4 

34TV 2 + 33TV - 74 2TV 2 — 37TV - 78 2 14TV 2 + 13TV - 34 2 

-O3O1 + —. ■ . !v\ 2 ~ A AT/ A T -1 \ / AT 2>->l 



3TV(TV + l)(TV + 2) 8TV(TV + l)(TV + 2) 1 4TV(TV + 1)(TV + 2) 

10TV 2 + 21TV + 6 TV 2 -TV-4 / x 

~2TV(TV + 1)(TV + 2) S2C2 + 8 (TV+l)2(TV + 2) V 2,1 + ^ 3 + 2S ~ 2Sl ) 

5 TV 5 + 10TV 4 - 20TV 3 - 62TV 2 - 40TV - 8 2TV(TV + 3) 

TV 2 (TV + l) 2 (TV + 2) 2 2,1 ~ 3(TV + 1) 2 (TV + 2) C3 

35TV 5 + 66TV 4 - 118TV 3 - 294TV 2 - 152TV - 40 

+ 3TV 2 (TV + l) 2 (TV + 2) 2 3 

TV 5 + 6 TV 4 + 4TV 3 - 30TV 2 - 40TV - 8 

(201 



TV 2 (TV + 1) 2 (TV + 2) 2 
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2iV 2 (iV + l) 2 (iV + 2) 2 
N 5 + QN 4 + AN 3 - 30N 2 - AON - 8 3 iV 4 + 2iV 3 + 7N 2 + 22N + 20 

S l ~ 8 /AT" , ^3/ AT" , 0\2 6 



6iV 2 (iV + l) 2 (iV + 2) 2 1 (iV + l) 3 (iV + 2) 

Ao^ 2 - AiS? 2iV 4 + lliV 3 + 15iV 2 + 12iV + 8 

"I /at- , , r,\9 C2 



-2 



27V(iV + l) 3 (iV + 2) 3 (iV + l) 3 (iV + 2) 2 

P22S1 P23 



(79) 



iV(iV + l) 4 (iV + 2) 4 (iV + l) 5 (iV + 2) 4 

6iV 6 + 36iV 5 + 60iV 4 - 99iV 3 - 390iV 2 - 316iV - 40 , 
2N 6 + 20iV 5 + 40iV 4 - 45iV 3 - 170iV 2 - lOOiV + 8 , 

AN 8 + 50N 7 + 22AN 6 + 5AAN 5 + 927 TV 4 + 11407V 3 + U2N 2 - QAN - 208 , 
8N 8 + 8QN 7 + 370iV 6 + 805iV 5 + 807iV 4 - 16iV 3 - 772iV 2 - 568iV - 96 . 

{I / 4 9 2 2 
; ( 85*2 1 1 — 8S-i 1 — 55*4 — &Si So ^ — —S-tS-i — — Sn — S0S1 
2N(N + 2) V ' ' ' ' 3 2 2 1 

-I s ' - ^ - ^ - |*<0 + wliK' + s ^ + I s ' + 2S ^) 



N 2 + 7N + 8 

rC3 



+ 



3(iV + l) 2 (iV + 2) 2 

1 14iV 3 + 41iV 2 + 51N + 36 12iV 5 + 124iV 4 + 472iV 3 + 817iV 2 + 641iV + 176 
3 N(N + 1) 2 (N + 2) 2 3 2N(N + l) 3 (iV + 2) 3 2 

4iV 4 + 16iV 3 - 4iV 2 - Q1N - 48 2 iV(lliV 3 + 56iV 2 + 92iV + 49) 

+ 2N(N + l) 2 (iV + 2) 3 1 + (N + l) 3 (iV + 2) 3 C2 

^ + ,„. ±,,„J . (80) 



iV(iV + l) 3 (iV + 2) 4 (iV + l) 5 (iV + 2) 5 
16iV 6 + lOliV 5 + 194iV 4 + 4iV 3 - 421iV 2 - 501iV - 192 , 

62iV 8 + 668iV 7 + 3073iV 6 + 7849iV 5 + 12052iV 4 + 11127iV 3 + 5640iV 2 + 1065iV - 128 

TpC a \ ; r~, r f 45*2 1 1 — AS^ 1 — 45*2 1 5*1 — —SuSi — —SoS^ — — St — S^C^ 

1 2N{N + l)(iV + 2) V 2,1,1 3,1 2,1 3 2 1 12 1 1S 

4 n \ iV 3 - 17iV 2 - 41iV - 16 / ^ 00 l n , 

UN + 20 
+ AN(N + 1)(N + 2) 4 



7iV + 36 2 3iV + 4 1 lliV 3 + 177V 2 — iV + 16 

"8iV(ivTT)(ivT2) 2 + 2iV(iV + l)(iV + 2) 2 ^ ~ 3 N(N + l) 2 (iV + 2) 2 3 



N + A 10iV 5 + 48iV 4 + 122iV 3 + 222iV 2 + 213iV + 64 

rC3 H „,,/», ; lW » r ; ^ ^2 



3(7V + l)(iV + 2) 2 ^ 2iV(iV+ l) 3 (iV + 2) s 

2iV 5 + 48iV 4 + l74iV 3 + 242iV 2 + 1617V + 64 2 4iV 3 + 267V 2 + 517V + 32 
+ 2/V(iV + l) 3 (;V + 2) 3 1 + (JV + l) 2 (iV + 2) 3 C2 



P2%S\ P; 
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iV(iV + l) 4 (iV + 2) 4 (iV + l) 4 (iV + 2) J 



(81) 
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P 26 = 8N 7 + 120iV 6 + 595iV 5 + 1538iV 4 + 2432iV 3 + 2419iV 2 + 1329iV + 256 , 
P 27 = 22N 7 + 266iV 6 + 1360iV 5 + 3826iV 4 + 6400iV 3 + 6376iV 2 + 3515iV + 832 
A Q9 T C f 2 3 ^ 3 + C3 l ^ + N 2 -3N-1 ( \ 



s 
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3iV 2 (iV + l) 2 N 3 (N + l) 3 

+ N 5 {N + l) 5 (iV + 2) 2 } ' ^ 
AN 9 + 8N 8 + QN 7 + 35iV 6 + 66iV 5 - QN 4 - 85iV 3 - 61iV 2 - 24iV - 4 . 



3(iV-l)iV(iV+l) 2 (iV + 2) 



35 3 + Cs) 



2Af4 + ^ - ^ - 20jV - 16 <2S 2 + b) - > ,(83) 



(N -1)N(N +1) 3 (N + 2) 2 V *V (AT — 1)N(N + l) 5 (iV + 2) 4 
P 29 = 2iV 8 + 10iV 7 + 28iV 6 + 91iV 5 + 213iV 4 + 160iV 3 - 248iV 2 - 512iV - 256 . (84) 

C2 ^ I rri\ }\ If om ( m 2 \ ( 1 iV 2 + 3iV + 2 

i=i 



(N + l)(N + 2) \2(N+l)(N + 2) 8N 



i=i 



X ^-^ +£ (i + fjJ (A f + l)(A r + 2) ^ (86) 

Note that for diagram u and v the sum runs over all heavy quark flavors. 
The diagrams contributing to the pure singlet contributions yield 



A Qg 



T F cJ- 4 - {N + 2){N - 1) hs 3 + C.) - 2 5N3 ~ 5N2 - 16N - 4 (2S 2 + C2) 
F F \ 3 N 2 (N + 1) 2 \° D3 ^^J z AT3(Ar + i)3(Ar + 2 ) \ 2 ^ ^ 

+ iy 5 (Ar + l) 5 (Ar + 2) 3 } ' (87) 
p 30 = 5iV 9 + 25iV 8 -35iV 7 - 229iV 6 -107A^ 5 + 481A^ 4 + 688A^ 3 + 360iV 2 + 112iV + 16 . 

4 Q9 Tr f 32 3 ^ 3 + C3 32 2iV + 3 (2? +r\ 

b " F F [ 3 (N-l)N(N + l)(N + 2) {N-l)N{N + l) 2 {N + 2) 2 \ 2 

(2N + 3)(iV 4 + 6iV 3 + 5iV 2 - 12N - 16) 

iV(iV- 1)(N+ l) 4 (iV + 2) 4 ' ' 
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The flavor non-singlet contributions are given by 

, 2 (iV + 2)(iV-l) _ 2 iV 4 + 2iV 3 -10iV 2 -5iV + 3 
! " M 9 N(N + 1) ^ 3 9 iV 2 (iV + l) 2 ^ 2 

P31 



81 iV 4 (7V + l) 4 J ' 

P 3 i = 49iV 8 + 196iV 7 - 83iV 6 - 533iV 5 - 374iV 4 - 59iV 3 - 33iV 2 + 9N + 27. 



(89) 



656 „ 392 1 

-^ Sl + sr) • (90) 
TS?» = r F c4-|-l^ (91) 
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B Infinite Sums 



In this appendix we list a series of infinite sums which were needed in the present analysis 
and are newly calculated. In addition we made use of the sums in [7]. o\ is a symbol for 
Z)fcli(V^) an d t ne corresponding sums are divergent. The calculation was partly performed 
using integral representations, solving difference equations and using the summation package 
Sigma, see also [48]. 



B.l Weighted Beta functions 



E 



B(N,i) 
^{i + N + 2f 



N 



AS^ 2 {N + 2) + 2S_ 3 (N + 2) + 2C 2 S 1 (N + 2) + 2( 3 



N(N + l)(N + 2) 



-GS- 2 (N + 2) - 3( 2 
N(N + 1)(N + 2) 



+ 



N(N + 1)(N + 2) 2 ' 



(92) 



E a 



B(N,i) 
^(i + N + 3y 

B(N,i) 



N 2S. 2 (N + 3) + ( 2 N 4 + 5N 3 + 13N 2 + 18N + 13 

"v / at/at , 1\/Ar , n\/sr , n\ ' at f at , iN9/at , o\9/Ar , o\9 ' 



iV(iV + l)(N + 2)(N + 3) N(N + 1) 2 (N + 2) 2 (N + 3) 5 



E 



= 24(-l) 



N 



2S_ 2 (N + 4) + ( 2 



+ 



N(N +1)(N + 2)(N + 3)(N + 4) 

N 6 + 11N 5 + 54iV 4 + 143iV 3 + 213iV 2 + 178iV + 100 
N(N + 1) 2 (N + 2) 2 (N + 3) 2 (N + A) 2 

-l + 2iV -1 + N + N 2 

- + (1 - N)S lt2 (N) + + AT + S 2 (N) 



(93) 



(94) 



+ 



N 

1-N-N 2 
TV 



N 



+ (-l + N)S 1 (N)y 2 + (l-N)( 3 



(95) 



E 



B(N,i) 
i (< + 4) 



-144 + 300iV - 415iV 2 + 241iV 3 - 63iV 4 + 6iV 5 
144AT2 



B.2 Weighted Beta functions and harmonic sums 



(96) 



E-rr 5iW 

i=i 



(N - 1)S 3 (N - 1) - (N - 1)C 3 - S 2 (iV - 1) + C 2 , (97) 

N 2 - 3iV + 3 



+ 



2 -TV 



(98) 



= (-C 3 + 5 3 (iv-i)) 



(iV — l)(iV — 2)(iV — 3) 
6 
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i=i 



i + 4 



i=l 



i + N + 3 



E ? 



B(AT,i) 



1=1 



i + N + 4 



Si(i) 



E 



B(N,i) 
^(i + N + 1) 



Mi) 



1=1 



f£(M 5l(j + JV ) 



j=i 



i + 2 



/ 1 ^ 3iV 3 -l8iV 2 + 33iV-22 (6iV - 13)(iV - 3) 
l C2 " Ss(iV " 1} J 12 + 24 



(99) 



+ 



-864 + 2040iV 2 - 3746iV 3 + 2453N 4 - 675iV 5 + 66iV 6 

864/V3 

300 - 550iV + 385N 2 - llOiV 3 + 1UV 4 



144 

(N - 1)(N - 2)(N - 3)(N - 4) 



((2-S 2 (N) 
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(Cs - S 3 (N)) 



(100) 



C 2 - S 2 (N + 2) 3N 6 + 15N 5 + 36N 4 + 51iV 3 + 52N 2 + 3QN + 8 



N + 3 



N 3 (N + l) 3 (N + 2) 3 



(101) 

C 2 - 5 2 (iV + 3) 2N 9 + 2AN 8 + 129iV 7 + 408N 6 + 854iV 5 + 1270N 4 



+2 



N + 4 " N 3 (N + 1) 3 (N + 2) 3 (N + 3) 3 

1405iV 3 + 1158iV 2 + 594iV + 108 
N 3 {N + 1) 3 {N + 2) 3 {N + 3) 3 ' 



-1) 



N 



2S 1 ,- 2 (N)+S- 3 (N) + ( 2 S 1 (N) + ( 3 ( 2 -S 2 {N) 



N(N + 1) 



+ 



(N+lf 



(N - 1) (2S 3 (N) + S 1 (N)S 2 (N) - (^(N) - 2C 3 ) 

1-N + N 2 



N- 1 



(5 2 (JV) - C2) 



+ 



AT2 



(102) 



(103) 



(104) 



(N 1} 2 (iV 2) (-25 s (iV) - S 1 (N)S 2 (N) + C 2 ^(iV) + 2C 3 ) 



i=i 



(N-2)(2N-1) 
2N 

2-2N 2 + N 3 



2N 3 

(N — l)(N — 2)(N — 3) 
6 



(S 2 (N) - C 2 ) - 



-4 + - 7N 2 + 2N 3 
AN 2 



^i(iV) 



(105) 



2S 3 (N) + S 1 (N)S 2 (N) - C 2 Si(N) - 2£ 



i=i 



_ (j V-3)(3^ 6j V + 2) (s2(jv) _ C2 

-12 + 21N 2 - 19N 3 + AN 4 
12A^ 

36 - QQN + 85 N 2 - 39 N 3 + QN 4 n . _ 

+ 36iV^ SliN) 

144 - 340N 2 + 397N 3 - 150N 4 + 18N 5 
144AT3 

-144 + 300N - 415N 2 + 241N 3 - 63N 4 + 6 A" 5 
UAN 2 



(106) 



Si (AO 
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1=1 v ' 



2 r (i + jV + 3) g .(' + Af ) 



^B(iV,i)S' 1 (i + AO 2 = 



i=l 



f;^M Sl(i + JV) » 

i=i 1 



+ 
+ 



(N - 4){2N - 3)(1 - 3JV + N 2 ) 
127V 

(JV - 1)(JV - 2 ^- 3)(JV - 4) gl w} ( 6 - w) 

(j V_ 1)(J v_ 2)(Af _3) (j y-4)^_ g3( ^ 



12 



(107) 



y^(45 1> _ 2 (N + l)-25_ 3 (iV + l) 



iV(iV+ 1) 

-25x(iV + l)S_ 2 (iV + 1) + C 2 ^(iV + 1) - Cs - 2S_ 2 (iV + 1) 



- 6 ) 



^(TV + l) 



1 



+ N(N + l) 2 + iV(iV + l) 3 ' 



(108) 



-1) 



7V+1 



&S!- 2 (N + 2) - 4S_ 3 (iV + 2) 



iV(iV+ l)(JV + 2) 

-45! (JV + 2)S_ 2 (iV + 2) + 2C 2 5x(iV + 2) - 2( 3 

\ 1 + N + N 2 

- 10S _ 2(N+ 2)-^) + mN + +1) ; (N+2) 

1 + 7N + 6N 2 + N 3 
~ N(N + 1) 3 (N + 2) 3 ' 
16iV 3 + 12 + 40N + 30iV 2 + 6N 4 + N 5 



2 Si(iV + 2) 



(109) 



+ 



N 3 (N + 1) 2 (N + 2) 2 (N + 3) 

85iV 2 + 36 + 66iV + 69iV 3 + 34iV 4 + 9N 5 + N 6 
N 2 (N + 1) 2 (N + 2) 2 {N + 3) 2 

v( C(2) + 2S_ 2 (A0 



S'i(JV) 



-6(-l) 



(110) 



N(N + 1)(N + 2)(N + 3) 

144 + 564iV + 564iV 2 + 361N 3 + 180N 4 + 62iV 5 + 12iV 6 + N 7 
N 3 (N + 1) 2 (N + 2) 2 (iV + 3) 2 (N + 4) 

424iV 5 + 110N 6 + 16N 7 + N 8 



N 2 (N + 1) 2 (N + 2) 2 (N + 3) 2 (N + A) 2 



Si(N) 



+ 



576 + 1200N + 1660N 2 + 1576N 3 + 1013A^ 
N 2 (N + 1) 2 (N + 2) 2 (N + 3) 2 (N + 4) 2 

C(2) + 2S_ 2 (JV) 



Si(N) 



-24(-l) 



N : 



N(N + 1)(N + 2)(N + 3)(N + 4) 



(111) 



—2 + 5A" — 2iV 2 + iV 3 (1 - N + N^S^N) S^N) 2 
(N- 1) 3 N 3 + 2 (N - 1) 2 N 2 + N- 1 

(N-1)N (N — 1)N y ' 

^ + 2^ + _ 5 2 (iV)5 2 (iV) - 4S l( A0S 3 (A0 



-3^ 4 (iV)-2 ( - 1) ^- 2(iV) + 2^ 2 (iV) 



-1)% 
N 2 
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i + N + 1 



i=i 



i + iV + 2 



E 



i + N 



B(JV,i) 



1=1 



i + iV + 1 



i=i 



i + N + 2 



B(JV,i) 



i=l 



i + N + 2 



+S 2 (N)( 2 + 25_ 2 (iV)C 2 + ^C 2 + 45i(7V)C 3 , 

l)^ 1 /, ^l(^) 2 ^(AH 2 

<( 2 + 2S- 2 (N)) + L +2^^ + 



AT2 

(-1) 



iV 2 



iV 3 iV 4 



(113) 
(114) 



N(N + 1) 



(-6S- 2 ,i(N - 1) + 2S_ 2 (iV - l)5'i(iV - 1) 
+S_ 3 (iV - 1) + C 2 5i(iV - 1) - 3C 3 - 2S_ 2 (iV - 1) - C 2 ) 



AT 2 (iV + l) 2 

2 + 37V 
"iV 3 (iV+ l) 2 ' 

(-1) N 



iV 2 (iV + 1)< 



(115) 

-12S'_ 2i i(iV - 1) + 4S_ 2 (iV - 1)51 (iV - 1) 



N(N + 1)(N + 2) 

+2S„ 3 (N - 1) + 2C 2 5i(iV - 1) - 6Cs + 2(N - 1)S- 2 (N - 1) 

, sr c _\ 4 + 6iV + 7iV 2 + 4iV 3 + iV 4 /Ar _ 
+(7V - 5)C 2 j + Ar2 , Ar , TW7? , n , 2 5i(iV - 1)' 



4 + 6iV + 7iY 2 + 4iV 3 + N 4 
N 2 (N + 1) 2 (N + 2)2 l( 



7V 2 (7V + l)(7V + 2) 5 



7V 3 (7V + l)(iV + 2)2 ' 



Si« 2 = 



5i« 2 = 



g^CjV - 1) - 25 3 (iV ~ 1) ~ C 2 £i(iV - 1) + 3C 3 

N 

S lt2 (N)-2S 3 (N)-( 2 S 1 (N) + 3( 3 S 2 (N)-( 2 , 2 



(116) 
(117) 



iV+ 1 



N 2 



+ iV 4 ' 



(118) 



Si, 2 (N + 1) - 2S 3 (iV + 1) ~ C 2 £i(iV + 1) + 3C 3 
N + 2 



S 2 (N + 1)-C 2 ) 



+ 



2iV 2 + + 1 
' iV 2 (jV + l) 2 

5iV 4 + 8N 3 + 13iV 2 + 8iV + 2 



S 1 (i)S 1 (N + i) = 



iV 4 (iV + l) 4 



(119) 



iV(iV + l)(JV + 2) 



45_ 2jl (iV) - 65_ 3 (AT) - AS^^S^N) 



-2C 2 5x(iV)-2C 3 -2 



(2 



-4 



- 2 



5 3 (iV) 
'N + 2 



S 1 (N)S 2 (N) + ( 2 S 1 (N)_ + 2C 3 



iV + 2 iV + 2 iV + 2 

2 + 7iV + 7jV 2 + 5iV 3 + ^4 

+ N 3 (N + l) 3 (iV + 2) l( } 
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^ B(N, i)S 2 (N + i) 



i=l 



E 



i=i 



N + i 



i=i 



iV + i + 1 



i=l 



iV + i + 2 



i=i 



iV + i + 2 



+2 



2 + 7N + 9N 2 + 4N 3 + N 4 
N*(N +l) 3 (N + 2) ' 

1 , S 2 (N) 2 (-l) N S- 2 (N) 



(120) 



-1)% 



(N-1)N 2 N-l 



(N-1)N (N-1)N 



(121) 

^^-S 2 (N + i) = -3S A (N)-2S„ 2 (N) 2 -S 2 (N) 2 -2S„ 2 (N)( 2 + S 2 (NK 2 + 1-(I 
i 10 

(122) 
(123) 



+ (-1) 



A? 



C(2) + 2,S_ 2 (iV) , S 2 (N) 



N 2 



+ 



N 2 



-1) 



A' 



C(2) + 2,S_ 2 (iV) , S 2 (N) 



N 2 



+ 



N 2 



N 2 (N + iy 

S-s(N) 



2S lt . 2 (N) 2S^ 2 (N) 



f4 



N(N + 1) 
Si(N)( 2 Cs 



N 2 



N(N + 1) N 2 N(N + 1) N(N + 1) 
N 4 + AN 3 + 7N 2 + 6iV + 4 



(124) 



+ (-1) 



iV 2 (iV + l) 2 (iV + 2) 2 
45 1> _ 2 (iV) 



-S 2 (N) + 



N(N + 1)(N + 2) 
25_ 3 (iV) 



iV(iV + l)(iV + 2) 

iV(iV + l)(iV + 2) 
5 1)2 (iV) 5x(iV)C 2 



+ c 2 



iV(iV + l) 2 (iV + 2) 

2(iV - 2)S- 2 (N) 
N 2 (N + 2) 

2S 1 (N) 



N-2 



N(N + 1)(N + 2) iV 2 (7V + 2) 



(125) 



+ 



Ca 



(2 + 3iV + 4iV 2 + iV 3 ) 1 S 2 (iV) 



iV + 2 N + 2 N + 2 
(2 + 3iV + 4iV 2 + iV 3 )C 2 



iV 2 (iV + l) 2 (7V + 2) 
2 



N 2 (N + 1) 2 (N + 2) N(N + 1) 4 (N + 2) 



(126) 



-8 + N + 3N 2 



S 2 (N) 



+ 



C-2 



N(N + 1) 4 (N + 2) 2 (N + 2) 2 (N + 2) 2 

2(2 2S 1 (N)( 2 

N(l + N) 2 (N + 2) N(N + 1)(N + 2) 

2Cs 4S_ 2 (iV) 



+ 
+ 



+ 



N(N + 1)(N + 2) N(N + 1) 2 (N + 2) 
2S- 3 (N) 4S , 1 ,_ 2 (JV) 



+ 



N(N + 1)(N + 2) N(N + 1)(N + 2) 



)(-D 



A? 



(127) 
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B.3 Weighted Harmonic Sums 

~' S 1 (i)S 1 (i + N) 



i + N 



E 

i=i 

^S } (i)S 1 (i + N) 



i=i 

oo 

E 

i=i 



E 



E 

i=i 



i + iV + 1 

S 1 (i)S 1 (i + N) 
i + N + 2 



S 1 (i)S 1 (i + N) 
i + 3 



i + 4 



E 

i=i 



E 

i=i 



S^i + N)^) 



S 2 1 (i + N)S 1 (t) 



E 

i=l 



E 

i=i 



S^Szii + N) 



S 1 (i + N)S 2 (i) 



E 

i=l 



i + N 



d ~ & - *?f> - 2 ^ + ± ^ - gl (JV^JV) , (128) 

(j?-4C3 , 5 3 (7V) - 5?(JV) S?(iV) 



(129) 



+ 



7V + 1 



Si(A% 



c 2 



7V + 1 



(iV + 1) 2 



(iV + 1) 2 ' 

3iV 2 - 6iV + 2 



(130) 



3 ' ' N(N — 1)(N — 2) 

(TV + l)(7iV — 6) 



5i,i(7V) + 



13JV-19 



-Si(JV) , 



4N(N - 1) 

1 3 (N+ l)(132 - 232iV + 85iV 2 
3^1 - - 



4(iV- l)(7V-2) 
(131) 



36(iV - 2)(JV - l)iV Sl{N) ~ 3 Cs + 5l ' 1 ' l(7V) 



2 (-3 + 2iV)(l-3iV + iV 2 ) ^ + 809 - 909iV + 232iV 2 



(N -3)(N -2)(N - 1)N 



36(N - 3)(N - 2)(N - 1) 

(132) 



^ + |c 2 2 + 3^(iVK3 + 

, S*(N)S 2 (N) , 2 



S*(N)-S 2 (N) 



C 2 -5'i(iV)52 1 i(JV) 



(133) 



of 43 
4 20 



C 2 2 + 5S 1 (A0C3 + 



35 2 (7V) - S 2 (iV) 



C2-2S 1 (N)S 2 , 1 (N) 



+S 2 1 (N)S 2 {N) + S^SziN) - 



Si(N) , Sf(N) 



+ 



(134) 



S*(N) - 3S 2 (N) 



2 C 2 -25 3il (iV) + -5 4 (iV) 



Sf(N)S 2 (N) 



+ S 1 (N)S 3 (N) , 



(135) 



= f C2 -f + ^)C3 + ^^ 



C2 + 25 , i,i )2 (iV) 



-25 1)3 (N) + S 1 (N)S 2 , 1 (N) - Sf(N)S 2 (N) + S ^ N) - S ^ N) , (136) 



at a 2 



i + E1 C2 _ _L C | _ + ?M1 _ ?M1 _ 3 ^{N)S 2 (N) 



9 



8 4 ' 40 



+ 



2N 



-3 



8 



+ *(*)(3^-*<lV)) + ^, 



(137) 
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^ jUi+jWj) ^ , , 39 2 S*(N)+S 2 (N) 



=1 



8 4 40 

+3£i (^ + S^V) _ S|^V) + Si(jV)( _ S2i(Ar) + sm) 
^ ■ (138) 

S 2 (iV) S 2A (N) , 0/A7 .f n S 2 (N) 



2 JV 



+ Si(JV) (2^^ + 5 2;1 (7V) - 2S 3 (iV)) 



+ 2^ + ^V)-?*f>, (139) 



2- — TTiv — - Y C2- To C2 + liv" 25l(iV V C3+ lnv 2 2~ ) (2 

i=l 



s K n) s 2{ n) + sm)_s^n) + Si(n)SjAn) + SsAN) 



N 2 N 2 2 N 

SAN) , x 

+^Y 1 , (140) 



^^i^+AO^) 6 2 / 2 x ^(JV) 5 2 (N) SI(N) ^f) 
^ (i + N) 2 ~ 5 C2 + [~N +2Sl{N) ) C3 + ^ 2 ~ + ^ 2 2~ + ^\T~ 



-5i(JV)5 2>1 (JV) - 5 3>1 (7V) - , (141) 



^ S 1 (i + N)SHi) o\ 3( 2 ( 2 oq(m \ r ^fSi(N) S 2 (N) S 2 (N)^ 

4- ^+iv " T"X + liv" 25l(iV V C3+ l^v 2 2~ J C2 



+ S 3 (iV) Sf(JV) , c2 , Ar ^ 1 35 2 (iV) A S 2 (7V) S 2 (iV) 



AT 4 V iV 2 2 J N 2 



v + (3^U + 5 2>1 (JV) - 25 3 (N)) + + 5 3>1 (7V) 

-5 4 (iV), " ' (142) 



B.4 Harmonic Sums 



^S 3 {i + N)-S 3 (i) = S 2 (N)-(N + 1)S 3 (N)+N( 3 , (143) 
i=i 

f; Wj + J v)-s,(j))i = s.W-s.w + ^ + i)s,(Af)-^ + i) 6 



i=i 



-NS 2 (N) + NC 2 , (144) 
^(s^i + JV)-^*)) 3 = -ls*(N)-S*(N)-^S 2 (N) + 3NS 2tl (N)-NS 3 (N) 

+N( 3 , (145) 
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i=i 



J2(s 1 (i + N)-S 1 (i)yi = ?L(-(N + 1)( 3 -3(N+1)S 2A (N) + (N + 1)S 3 (N) 

i=l 



H 3N-l )C2 ) + 1 + 2 \ +6N2 S 2 (N ) + lsm 

+ 2 -^ Sl (N) , (146) 



J2( S ^ + N ) -Si(i))s 2 (i) = --S 2 1 (N)--S 2 (N)+NS 2 , 1 (N) + N( 2 -NS 1 (N)( 2 

+Na 1 ( 2 - 2N( 3 , (147) 
^(s^i + ^-S^Szii + N) = -(1 + S 1 (N))S 2 (N)+NS 3 (N)-NS 1 (N)C 2 

+N(( 2 + a l ( 2 -( 3 ) , (148) 



i=i 



B.5 Miscellaneous Sums 



fc=0 



iV - 1 /AT - 1\ 

£ , )(-l)'S(i + l-e/2,2- e /2) = S(l-e/2,JV + l-e/2). 

B.6 Double and Other Sums 

^ S 1 (i)S 1 (i+j + N) S 1 (N) { SfjN) g 2 (iV)x g(AQ S*(N)S 2 (N) 

4^ i(i + ?)(? + AT) "~ iV ^^^l jy N ) QN N 



(150) 



^ i(i+j)(j + N) 



N N 



+ , l( ,)(_3« + 4^)) 



Ai(AQ S 4 (1V) 
~ 2 ^V 2iV~ ' U5ij 



^ N + k ~ iV2 
fc=i 



g J yv 2g_ 2 (iV) + C 2 gi(AQ 
2[ l ' N N 2 

J ( u n S-2,i(N) n 2S^ 3 (N)-( 3 
+ \ { } N + 1 ' AN 

+ 4iV Sl{N) + J 
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+ 



48N 



y ("1)% 2 - 6(~l) N S 1 (N)( 3 + ( 2 (3(-l) N S^N) 



E 

fc=i 



B(iV + l,fc + e/2) 
iV + fc 



+3(-l) JV S 2 (JV)) + 48(-l) JV S_ 2 , 1 , 1 (JV) - 24(-l) JV 5_ 3 ,i(JV) 
-24(-l)^_ 2 , 2 (A0 + ^(-l)^^) + 6(-l) JV ,S 1 2 (iV)^ 2 (iV) 

+6(-l) JV S 2 (JV)S- 2 (JV) - 2- 5s(JV) ^ 



+S' 1 (iV)(-3 
+0(e 4 ) . 



S 2 (AQ 
N 



N 



2A(-1) N S_ 2A (N) + 12(-1) N S- 3 (N) 



N 



N , 



TV 



i)» J2(- 1 Y ( E — fc + £/2) + 1 + e/2) 



,fc=i 



|( 1)jy ^ B(l,fc + e /2) ^ 



fe=i 



(152) 



(153) 



E 

fc=i 



£(A; + e/2,iV+l) 
N + k 



If 2S_ 2 (N) + ( 2 



+ -C 3 + C 2 5i(iV) + 25 li _ 2 (JV)-25_ 2)1 (JV) 



-C 2 2 -C 3 ^i(iV) + C25i,i(iV) 



+2{s , lil ,_ 2 (iV) + S- 2 ,i,i(N) - S lt - 2 ,i(N)} 



C 5 ^i(iv) WAO 

— H ^^<,2 3 (>3 H 3 <,2 



20 



VwW + Si,i,i,-2(JV) - S- 2tl , ltl (N) - S 1 , 1 ,. 2 , 1 (N) 



+0(e 4 ) . 



(154) 



B.7 Expansion of harmonic sums for small argument 

One may expand nested harmonic sums into Taylor series w.r.t. the outer argument, using the 
corresponding differentiation rules [51]. In the present calculation we made use of the following 
relations. 



S 1 (e) = ( 2 e-( 3 e 2 + - ( 5 e 4 + 0(e 5 ) , 

o 

S 2 (e) = 2C 3 e-^ 2 V + 4C 5 5 3 -^ 2 V + 0(5 5 ), 
Ss(e) = Ule - 6C 5 £ 2 + ^tfe 3 - 15C 7 £ 4 + 0(e 5 ) , 



(155) 
(156) 
(157) 
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S 4 (e) = 4C 5 5-yC 2 V + 20C 7 e 3 -yC 2 V + O(e 5 ) , (158) 

5 2> i(e) = ^+(2C3C 2 -yC 5 )^ + 0(e 3 ), (159) 

53,i(e) = (3(2)6 + 0(e 2 ) , (160) 

52,i,i(e) = (yC 5 -2C 3 C 2 )e + 0(e 2 ) . (161) 

These relations can be obtained expanding the representation of the sums in terms of Mellin 
transforms of Nielsen integrals weighted by 1/(1 ± x). In case of the single harmonic sums the 
expansions result from Euler's ^-function and its derivatives. 

B.8 Sample Calculation for one of the Sum 

In the following we illustrate the calculation of sum (151) in using the Sigma package : 

^ 1 ™ Sl (i)S 1 (i+j + N) 

First, we treat the inner sum for N fixed, 



By Sigma's creative telescoping algorithm we compute the recurrence relation 

-(J + N + l)fF(j) + (j + 1) (3j 2 + 3Nj + 7j + 2N + 4) F(j + 1) 
-0" + 2) (3j 2 + 3Nj + llj + 4N + 10) F(j + 2) 

+(j + 2)(j + 3)(j + N + 3)F(j + 3) = A(j) + B(j)S 1 (j + N) 

(164) 

where 



(f + 3j + 2) (j 2 + (2N + 3)j + N 2 + 3N + 2) 2 
x (V + (9N + 26)j 4 + (lOiV 2 + 63iV + 86) f + (5iV 3 + 49iV 2 + 156iV + 137) f 

+ (N 4 + 12N 3 + 74N 2 + 163N + 106) j + 5iV 3 + 34iV 2 + 61iV + 32^ , (165) 

n 3j 3 + (4iV + 13)j 2 + (-N 2 + UN + 18) j - N 3 - 2N 2 + 7N + 8 

{J) " (j + l)(j + 2)(j + N + l)(j + N + 2) ■ [ ] 

Next, we apply Sigma's recurrence solver and obtain three linearly independent solutions of the 
homogeneous version of the recurrence: 

1 S 1 (j + N) - 3 S 1 (j)+jS 1 (j + N) + l 
j" j J 2 (N + 1) 
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and one solution of the recurrence itself: 

PU) = + Nf - + + (-^ + S2(J + *° + —) S,(J + N) S * iN) 



6j 2f V 2j 2j Nj+jJ 3j 

Si(N) 2 37V + 2 (j-2(N+l))S 1 (N) S 2 (N) S 2 (j + N) S 3 (j + N) 
2Nj + 2j ■> f(N + l) 2 2f{N + l) 2 2j 2 2f 3j 

yj Si(i+N) yj S^i+N) yj S^Sxji+N) yj S 2 (i+N) yj S^i+N) 2 
| Z-<i=l j2 A-ii=l j Z-n=l i |_ Z-/t=l j |_ A-ii=l j 

j j 2 j 2j 2j 

+Sl0) {WTW wn~~ 3 ) 3 • ( ' 

The function F(j) is given by 

1 ^(j + AQ -jgiC7) + J^iC7 + iVT) + 1 
= «ij + «2 + a 3 j2(jv + l) + (168) 

for some properly chosen constants a±, a 2 and a 3 which are free of j. Looking at the initial values 
for j = 1, 2, 3 of F(j) we can conclude that 

= 4*w*-i^ + B*w-2civW)* (jv) -^ w+26 ' (i69) 

02 = U- Sm >- 2 4m + S2{ N ) + 2 J^P^ \ (170, 



2\ zv ' (N + l) 

«s = ^(N+l)S 1 (Nr + S 1 (N)- 1 ^ T . (171) 

One obtains 

m = + W + M _« Sl(iV) + ^ + iV ) 3 ^(j + ao 2 



6j 2j 2 \2jN 2 2j J v ' Qj 2j 2 

+Si{j)2 f S 1 (j + N) _ J_\ + S^V) _ S 2 (j + N) _ S^N)_ + Ss{j + N) 



2j 2jNJ 2f 2f 3j 3j 

^2,l7VJ . 2^=1 i+AT . <Li=l — p — U + iV )Li=l — i — ^»=l (i+iV) 2 



J 2j 2j j 2 iV 



v-j Ai i+ivr v J ^2 ^ /a/' , i\a 1 or \ or 



,3 



2 j 2j V 2 j 2j ; j 

(172) 



Finally, we look at the indefinite nested sum 



^ = ( 173 ) 

.7=1 J 
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with 



lim S(N, a) = ±± S f Sli ^. 3+ A !" ) . (174) 

At this point we emphasize that the sum expression (173) with the derived sum representation 
of F(j) fits into the input class of Sigma. Hence, we can apply Sigma's machinery again and 
arrive for S(a, N) at the following sum representation 

q(n m SiW , S 1 (a + N)S 1 (Nf (S 2 {a) S 2 (N) (2 \ ( . 2 
S{a ' N) ~ "^iV~ + ON + {-4N APT + 2Nj Sl{N) 

+ ( S (a + N)( S2{N) M S2 > l(N) I 2(3 ^ 9 (N) S ^ a + N ^ 
+ { S * a + N > {-^F-^p)-^^ + w) Sl{N) 24N~ 

fsr^a Si(i+N)\ 2 

W _ S 2 (a + Nf (E,=i — i— J 3 / S^a + N) J_\ 

8N 8N 2N 11 j V 3iV 3N 2 J 

S(a?( SlW I Sl{a + N) S2{N K 1 ^ 4- 9 (a) ( S2{N) + 1 ^ 

S 3 {N) S 3 (a + N) SJJV) S 4 (a + N) EUj^W g^jjw 
3iV 2 3iV 2 4iV 4iV 3iV 2iV 

/frfr) gi(a + AQ 1 \^S 1 ( i ) 2 DUjff E°=i^ 
2iV 2iV 2iV 2 J ^ i + iV 6N 6N 

S 1 (a) _ Si(a + N) \ ^ S^ + iV) 
2iV 2iV J f-' i 2 



8iV 



+ 
+ 
+ 



1=1 

5i(a) Si(a + iV) 1 ^A^i + JV) 

, "Hv - + iV + iV 2 J ^ (i + iV) 2 
x ' i=i v ' 

spa Si(i+jV) „ a S^Siii+N) V a SifflSift+AQ V^a Si(i) 2 Si(i+AQ 

l^i=i (i+jy)3 _ 2^=1 g , Z^i=i (i+jy)2 2^=1 i+jy 
N " 2N N " " "' N 

( S 1 (a) _ S 1 (a + N) _ J_\ S 1 (i + iV) 2 _ 
\ 2N 2N N 2 ) 2s i N 

oV^a Si(i+AQ 3 v^a Si(i+N)S 2 (i) 



I f '1 = 1 l I f '»=! I I 

3.V 27V 

v^a 5i(i)5i(t+Af) 2 spa Si(i)S 2 (i+N) spa Si(i+N)S 2 (i+N) 
I Z-/j=l j Z-<i=l j |_ Z-u=l j 



gi(a) _ ^i(a + iV) \ A S 2 (i + N) 
2N 2N ) 2s i 

' i=l 



2N 2N N 

S 2 (a + N)(l-3N 2 ( 2 ) S 2 (N) (3N 2 ( 2 - 1) Si (a + N) 

6iV 3 + 6iV 3 + A^ 2 

[2s — i — ) n W 2 ) Sl{a) ~^m 

S 2 (a) S 2 (N) C2 1 \ , N)2 ( l-N 2 ( 2 S 2 (a + P 
^N~2PP +{) 



2N 2N ' N 2N 3 ) ' ' " v V 2A^ 3 4iV 
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+Si(« + AO ^— ^^ + ^^- aTJ 

^ " ~6i^ + 2iV + \2N* ~ ~W ) Sl{N) 

+5l(a + iY) ^^AT"2A^J ~^AT + ^iV iv" + (175) 

We remark that all the sums in this expression are algebraically independent, i.e., no relations 
occur that could cancel some of the involved sums. 

Finally, we send a to infinity in the last expression and note that the involved sum expressions 
can be simplified by the sum identities (128)-(142) and some additional identities of similar type. 
In the final expression divergences of the type o\ being contained in some of the terms vanish. 
We find the right hand side of (151). 
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C The first moment of the operator matrix element 



After the analytic continuation from the even values of N to N e C is performed one may 
consider the limit N —* 1. In this procedure the term (1 + (— 1) )/2 equals to 1. At 0(a 2 ) the 
terms oc T f Ca contain 1/z contributions in momentum fraction space and their first moment 
diverges. For the other contributions to the un-renormalized operator matrix element (after 
mass renormalization to 2-loop order), the first moment is related to the Abelian part of the 
transverse contribution to the gluon propagator Hy(p 2 , m 2 )| p 2 =0 , Figure 3, except the term oc Tp 



-CD-- 




Figure 3: Abelian part of the gluon self-energy due to heavy quarks. 

which results from wave function renormalization. This was shown in [6] up to the constant term 
in e. One obtains 



U v (p 2 ,m z ) = S £ a s T F U { y\p 2 ,m z ) + S'y s C F T F U ( y'(p 2 ,m') + 0{a 



r( 2 )^2 



0, 



with 



lim Uy\p 2 , m 2 ) 
lim Uy\p 2 , m 2 ) 



-A 



(1),N=1 

Qg 



n A Qg \Cf 



(176) 

(177) 
(178) 



Here we extend the relation to the linear terms in e. For the first moment the double pole 
contributions in e vanish in (177,178). From the corresponding QED-expressions 11^' given in 
[59] by asymptotic expansion of the photon propag ator {l/p^U^ip 2 , m 2 ) in m 2 /p 2 and adjusting 
the relative color factor for k = 2 to 1/4 = 1/ (CpC^), due to the transition from QED to QCD, 
the comparison can be performed up to the constant term in e. One obtains 



1 ~ 



lim — n^Qo 2 , m 



p 2 ^o p 
1 ~, 

p 2 ^o p 2 



-A 



(1),N=1 



lim — Uy\p 2 ,m 2 ) 



2T F Qg 
1 



m 



2\ e/2 



Ar,„ \C F 



ni 



2\ £ 



2T F C F ^ \fji< 
The latter term is easily obtained using MATAD [61]. 



8 e 

7T" + 0C2 

3e 3 



4 /31 
- + 15- — 
£ V 4 



(179) 
. (180) 
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